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Abstract 

We present a pedagogical review of neutrino physics. In the first lecture we describe the theo- 
CN ■ retical motivation for neutrino masses, and explain how neutrino flavor oscillation experiments can 

probe neutrino masses. In the second lecture we review the experimental data, and show that it 
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I/-) , is best explained if neutrinos are massive. In the third lecture we explain what are the theoretical 

. implications of the data, in particular, what are the challenges they impose on models of physics 
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beyond the SM. We give examples of theoretical models that cope with some of these challenges. 
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I. INTRODUCTION 



The success of the Standard Model (SM) can be seen as a proof that it is an effective low 
energy description of Nature. We are therefore interested in probing the more fundamental 
theory. One way to go is to search for new particles that can be produced in yet unreached 
energies. Another way to look for new physics is to search for indirect effects of heavy 
unknown particles. In this set of lectures we explain how neutrino physics is used to probe 
such indirect signals of physics beyond the SM. 

In the SM the neutrinos are exactly massless. This prediction, however, is rather specific 
to the SM. In almost all of the SM extensions the neutrinos are massive and they mix. 
The search for neutrino flavor oscillation, a phenomenon which is possible only for massive 
neutrinos, is a search for new physics beyond the SM. The recent experimental indications 
for neutrino oscillations are indirect evidences for new physics, most likely, at distances much 
shorter than the weak scale. 

In the first lecture the basic mechanisms for generating neutrino masses are described and 
the ingredients of the SM that ensure massless neutrinos are explained. Then, the neutrino 
oscillation formalism is developed. In the second lecture the current experimental situation 
is summarized. In particular, we describe the oscillation signals observed by solar neutrino 
experiments, atmospheric neutrino experiments and long baseline terrestrial neutrino exper- 
iments. Each of these results separately can be accounted for by a rather simple modification 
to the SM. Trying to accommodate all of them simultaneously, however, is not trivial. In 
the third lecture we explain what are the theoretical challenges in trying to combine all the 
experimental indications for neutrino masses, and give several examples of models that cope 
with some of these challenges. 

These lecture notes are aimed to provide an introduction to the topic of neutrino physics. 
They are not meant to be a review. Therefore, many details are not given and many 
references are omitted. There are many textbooks [1] and reviews [2-4] about neutrinos. 
There is also a lot of information about neutrinos on the web [5, 6]. All these sources 
provide more detailed discussions with complete set of references on the topics covered in 
these lectures. Moreover, they also cover many subjects that are not mentioned here. 

In preparing the lectures I used mainly the recent review by Gonzalez-Garcia and Nir [4] . 
This review is a very good starting point to anyone who wants to learn more about neutrino 
physics. 
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II. NEUTRINO MASSES 



A. Fermion masses 

In general, there are two possible mass terms for fermions: Dirac and Majorana mass 
terms. All fermions can have Dirac mass terms, but only neutral fermions can have Majorana 
mass terms. Indeed, all the massive fermions in the SM, the quarks and charged leptons, 
have Dirac mass terms. The neutrinos, however, while massless in the SM, can have both 
Dirac and Majorana mass terms. 

Dirac masses couple left and right handed fields 

m D ip L ip R + h.c, (2.1) 

where mc is the Dirac mass and ipL and tpR are left and right handed Weyl spinor fields, 
respectively. Note the following points regarding eq. (2.1): 

• Consider a theory with one or more exact U(l) symmetries. The charges of i^l and 
ipR under these symmetries must be opposite. In particular, the two fields can carry 
electric charge as long as Q(^l) = Q(^>r)- 

• Since ipi and ipR are different fields, there are four degrees of freedom with the same 
mass, mo- 

• When there are several fields with the same quantum numbers, we define the Dirac 
mass matrix, {m D )ij 

Wyfelife)] + h.c, (2.2) 

where runs from one to the number of left (right) handed fields with the same 
quantum numbers. In the SM, the fermion fields are present in three copies, and the 
Dirac mass matrices are 3x3 matrices. In general, however, mo does not have to be 
a square matrix. 

A Majorana mass couples a left handed or a right handed field to itself. Consider ipR, a 
SM singlet right handed field. Its Majorana mass term is 

m M i^i; R , i> c = C^ T , (2.3) 

where ttim is the Majorana mass and C is the charge conjugation matrix [7]. A similar 
expression holds for left handed fields. Note the following points regarding eq. (2.3): 

• Since only one Weyl fermion field is needed in order to generate a Majorana mass 
term, there are only two degrees of freedom that have the same mass, m M . 
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• When there are several neutral fields, tum is promoted to be a Majorana mass matrix 

(rnM)ij(Wk)i(ipR)3- (2-4) 

Here, % and j runs from one to the number of neutral fields. A Majorana mass matrix 
is symmetric. 

• The additive quantum numbers of ip R and ip R are the same. Thus, a fermion field can 
have a Majorana mass only if it is neutral under all unbroken local and global U(l) 
symmetries. In particular, fields that carry electric charges cannot acquire Majorana 
masses. 

• We usually work with theories where only local symmetries are imposed. In such 
theories global symmetries can only be accidental. A Majorana mass term for a fermion 
field ip breaks all the global U(l) symmetries under which the fermion is charged. 

B. Neutrino masses in the SM 

The SM is a renormalizable four dimensional quantum field theory where 

• The gauge group is SU{3) C x SU(2) L x U{l) Y . 

• There are three generations of fermions 

Q L (3, 2) +1/6 , U R (3, l) +2/3 , D R (3, l)_ 1/3 , 

L L (l,2)_ 1/2 , E R (1,1)^. (2.5) 

We use the (qc,QL)q Y notation, such that qc is the representation under SU(3)c, Ql 
is the representation under SU{2)l and qy is the U(l)y charge. 

• The vev of the scalar Higgs field, H(l, 2) +1 / 2 , leads to the Spontaneous Symmetry 
Breaking (SSB) 

SU(2) L xU(l) Y ^U(l) EM . (2.6) 

The SM has four accidental global symmetries: baryon number (B) and lepton flavor 
numbers (L e , and L T ). It is also convenient to define total lepton number as L = 
L e + + L T , which is also conserved in the SM. An accidental symmetry is a symmetry that 
is not imposed on the action. It is there only due to the field content of the theory and by the 
requirement of renormalizability. If we would not require the theory to be renormalizable, 
accidental symmetries would not be present. For example, in the SM dimension five and six 
operators break lepton and baryon numbers. 

The SM implies that the neutrinos are exactly massless. There are several ingredients 
that combine to ensure it: 
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• The SM does not include fields that are singlets under the gauge group, N R (1, 1) . 
This implies that there are no Dirac mass terms {H)vZv R } (We define H = ir 2 H*.) 

• There are no scalar triplets, A(l,3)i, in the SM. Therefore, Majorana mass terms of 
the form (A)z/£z/^ cannot be written. 

• The SM is renormalizable. This implies that no dimension five Majorana mass terms 
of the form (H)(H)v l vl are possible. 

• U(1)b-l is an accidental non-anomalous global symmetry of the SM. Thus, quantum 
corrections cannot generate Majorana (H){H)v l vl mass terms. 2 

Both the neutrinos and the SU(3) C xU(l) em gauge bosons are massless in the SM. There 
is, however, a fundamental difference between these two cases. The symmetry that protects 
the neutrinos from acquiring masses is lepton number, which is an accidental symmetry. 
Namely, one does not impose it on the SM. In contrast, the photon and gluon are massless 
due to local gauge invariance, which is a symmetry that one imposes on the theory. This 
difference is significant. An imposed symmetry is exact also when one considers possible 
new physics at very short distances. Accidental symmetries, however, are likely to be broken 
by new heavy fields. 

One can understand the result that the neutrinos are massless in the SM in simple terms 
as follows. For a massive particle one can always find a reference frame where the particle 
is left handed and another reference frame where it is right handed. Thus, in order to have 
massive neutrinos, the SM left handed neutrino fields should couple to right handed fields. 
There are two options for such couplings. First, the SM left handed neutrinos can couple 
to right handed SM singlet fermions. This is not allowed in the SM just because there are 
no such fields in the SM. The second option is to have couplings between the left handed 
neutrinos and the right handed anti-neutrinos. Such couplings break lepton number. The 
fact that lepton number is an accidental symmetry of the SM forbids also this possibility. 

1 Our notation is such that we use capital letters to denote fields before SSB and lowercase Greek or Roman 
letters to denote the fields after SSB. While such differentiation is meaningful only to field that are charged 
under the SM gauge group, we use this convention also for SM singlets. 

2 Note that U(1)b-l has gravitational anomalies. Namely, once gravity effects are included the symmetry 
is broken and neutrino masses can be generated. The SM, however, does not include gravity. What 
we learn is that it is likely that the neutrinos are massive in any SM extension that include gravity, in 
particular, in the real world. 
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C. Neutrino masses beyond the SM 



Once the SM is embedded into a more fundamental theory, some of its properties men- 
tioned above are lost. Therefore, in many SM extensions the neutrinos are massive. Neutrino 
masses can be generated by adding light or heavy fields to the SM. (By light we refer to 
fields that have weak scale or smaller masses, while heavy means much above the weak 
scale.) When light fields are added the resulting neutrino masses are generally too large and 
some mechanism is needed in order to suppress them. New heavy fields, on the contrary, 
generate very small masses to the neutrinos. We are now going to elaborate on these two 
possibilities. 

There are two kinds of light fields that can be used to generate neutrino masses. First, 
right handed neutrino fields, N R (1, 1) , couple to the SM left handed lepton fields via the 
Yukawa couplings and generate Dirac masses for the neutrinos 

y N HL L N R + h.c. =>• m D WZu R , (2.7) 

where yjv is a dimensionless Yukawa coupling and the symbol indicates electroweak 

SSB. This mechanism offers no explanation why the neutrinos are so light and why the right 
handed neutrinos do not acquire large Majorana masses. Another way to generate neutrino 
masses is to add to the SM a scalar triplet, A(l, 3)i. Assuming that the neutral component 
of this triplet acquires a vev, a Majorana mass term is generated 

\ N ALIL L =>• m N 17[v L) (2.8) 

where \n is a dimensionless coupling. In that case it remains to be explained why the triplet 
vev is much smaller than the SM Higgs doublet vev. Note that the triplet vev is required 
to be small not only to ensure light neutrinos but also from electroweak precision data, in 
particular, from the measurement of the p parameter [8] . 

The other way to generate neutrino masses is to add new heavy fields to the SM. In that 
case the low energy theory is the SM, but the full high energy theory includes new fields. 
Using an effective field theory approach, the effects of these heavy fields are described by 
adding non-renormalizable terms to the SM action. All such terms are suppressed by powers 
of the small parameter v/M. Here v is the SM Higgs vev, which characterizes the weak scale, 
and M is the unknown scale of the new physics, which can be, for example, the Planck scale 
or the GUT scale. Neutrino masses are generated by the following dimension five operator 

^HHU[L L =► m„ = A^, (2.9) 

where A is a dimensionless coupling. 

Note the following points regarding eq. (2.9): 
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FIG. 1: The seesaw mechanism acquire its name from analogy with a children seesaw. The lower 
one child takes its arm of the seesaw, the higher the other child goes. For neutrinos, the larger Mn 
is, the lighter the neutrino becomes. (I Thank Oleg Khasanov for the figure.) 

• Taking the SM to be an effective field theory implies m v ^ 0. Since we know that new 
physics must exist at, or below, the Planck scale, it is rather likely that the neutrinos 
are massive. 

• m v is small since it arises from non-renormalizable terms. 

• Neutrino masses probe the high energy physics. 

• The neutrino acquires Majorana mass. Thus, total lepton number is broken by two 
units. 

• The situation can be generalized to the case of several generations. Then, genericaliy, 
both total lepton number and family lepton numbers are broken and lepton mixing 
and CP violation are expected. 

A famous example of a realization of this effective field theory description is the seesaw 
mechanism, see fig. 1. (Historically, it was invented within 5*0(10) GUT theory, and later 
implemented in many other models.) Consider one generation SM with an additional fermion 
singlet N R (1, 1) . The following two terms in the Lagrangian are relevant for neutrino masses 



where Mn 3> v is the Majorana mass of the right handed neutrino. After electroweak 
symmetry breaking the second term leads to a Dirac mass of the neutrino. In the {yi, v c R ) 
basis the neutrino mass matrix is 



-M N N C R N R + Y V HL L N R + h.c, 



(2.10) 





7 



where mo = Y v v. Using 



Tr(m„) = M N , \ det(m v )\ = m 2 D , (2.12) 
to first order in v/M N we find 

m VR = M N m VL = ^L. (2.13) 

M N 

Comparing this result with eq. (2.9) one identifies the new physics scale M with M N and the 
new physics coupling A with Y 2 . We learn that the seesaw mechanism is indeed a realization 
of the effective field theory approach for neutrino masses. 



D. Neutrino mixing 

Massive neutrinos generally mix. The neutrino mass terms break the accidental family 
lepton number symmetries. (Total lepton number is also violated if the neutrinos have 
Majorana masses.) This phenomenon is very similar to quark mixing in the SM. It is 
therefore instructive to describe both lepton and quark mixing in parallel. 

For quarks, the Cabibbo-Kobayashi-Maskawa (CKM) matrix, V, corresponds to non- 
diagonal charged current interactions between quark mass eigenstates 

fyui;)iVi j r(d L ) j W+ , i = u,c,t, j = d,s,b. (2.14) 

For leptons, it is common to use two different bases. The flavor basis is defined to be the 
one where the charged lepton mass matrix and the W interactions are diagonal. In the 
mass basis both the charged lepton and the neutrino mass matrices are diagonal, but the W 
interaction is not. In that basis the leptonic mixing matrix U is the analogue of the CKM 
matrix. 3 Namely, it shows up in the charged current interactions 

^hUecfMiW- , £ = e,n,r, i = 1,2,3. (2.15) 

When working in the mass basis, the formalism of quark and lepton flavor mixings are 
very similar. The difference between these two phenomena arises due to the way neutrino 

3 Recently, in many papers the matrix U is called the Maki-Nakagawa-Sakata (MNS) matrix, J7mnSj or the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, E/pmnSj since these authors were the first to discuss 
leptonic flavor mixing. Others call it the CKM or KM analogue for the lepton sector since Kobayashi and 
Maskawa were the first to discuss CP violation from such matrices. Here we adopt the notation of [4] 
and call U the leptonic mixing matrix. The fact that the CKM matrix is denoted by V helps in avoiding 
confusion. 
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experiments are done. While quarks and charged leptons are identified as mass eigenstates, 
neutrinos are identified as flavor eigenstates. Indeed, there are two equivalent ways to 
think about fermion mixing. For quarks, mixing is best understood as the fact that the W 
interaction is not diagonal in the mass basis. For leptons, we usually refer to the rotation 
between the neutrino mass and flavor bases as neutrino mixing. 

Note that the indices of the two matrices are of different order. In V the first index 
corresponds to the up type component of the doublet and the second one to the down type. 
In U it is the other way around. 

In general, U is described by three mixing angles and one Dirac CP violating phase. If the 
neutrinos are Majorana particles there are in addition two more CP violating phases called 
Majorana phases. The Majorana phases, however, do not affect neutrino flavor oscillations, 
and we do not discuss them any further. 

SM singlet states can mix with the three SM neutrinos. When the masses of the singlet 
states are small they can participate in neutrino oscillation. Then, the singlet states are usu- 
ally called sterile neutrinos while the standard model neutrinos are called active neutrinos. 4 
The mixings between all neutrinos are taken into account by enlarging the mixing matrix U 
to a 3 x (n + 3) matrix, where n is the number of sterile states. 

When the SM singlet states are very heavy (with mass M >> v) , they cannot participate 
in neutrino oscillations. Their couplings to the light neutrino is important since it gives 
them their masses. Thus, the heavy states are usually referred to as right handed neutrinos. 
The mixings relevant for neutrino oscillations is described by the matrix U that refers only 
to the 3x3 sub-matrix that corresponds to the mixing between the mostly active states. 
The presence of the heavy states affects the theory since in that case the matrix U does 
not describe the full mixing and, therefore, it is not unitary. This deviation from unitarity, 
however, is very small, of order v/M, and it is usually neglected. 

III. METHODS FOR PROBING NEUTRINO MASSES 

The conclusion of the last section is that it is very likely that the neutrinos have very small 
masses and that they mix. There are several ways to probe neutrino masses and mixing 
angles. Discussions about astrophysical and cosmological probes of neutrino masses can be 
found, for example, in [1, 3]. Here we briefly discuss kinematic tests for neutrino masses 

4 Oscillation occurs only between left handed states, that is, both active and sterile states are left handed. 
One should not get confused by the fact that the sterile states are left handed; the identification of left 
handed states as SM doublets and right handed states as SM singlets does not have to be true in theories 
that extend the SM. 
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and neutrinoless double beta decay. Then we develop the formalism of neutrino oscillation, 
which is the most promising way to probe the neutrino masses and mixing angles. 



A. Kinematic tests 

In decays that produce neutrinos the decay spectra are sensitive to neutrino masses. For 
example, in n — > fiu the muon momentum is fixed (up to tiny width effects), and it is 
determined only by the masses of the pion, the muon and the neutrino. To first order in 
ml/ml, the muon momentum in the pion rest frame is given by 

\P = 7, \ m 7T - m u 5 9 m v ■ (3-1) 

2m w \ n 11 ml - ml v ) 

Since the correction to the massless neutrino limit is proportional to ml, the kinematic tests 
are not very sensitive to small neutrino masses. The current best bounds obtained using 
kinematic tests are [8] 

m u < 18.2 MeV from r — > 5n + v , 
m u < 190 KeV from ir — > /iv , 

m u <2.2eV from 3 H ^ 3 He + e + V . (3.2) 

The sensitivities of neutrino oscillation experiments are much better than these bounds. Note 
that while oscillation experiments are sensitive to the neutrino mass-squared differences (see 
below) the kinematic tests are sensitive to the neutrino masses themselves. 



B. Neutrinoless double [3 decay 

Neutrino Majorana masses violate lepton number by two units. Therefore, if neutrinos 
have Majorana masses we expect there there are also AL = 2 processes. The smallness of 
the neutrino masses indicates that such processes have very small rates. Therefore, the only 
practical way to look for AL = 2 processes is in places where the lepton number conserving 
ones are forbidden or highly suppressed. Neutrinoless double (5 decay where the single beta 
decay is forbidden is such a process. An example for such processes is 

3 ^Ge -> % Se + 2e~ . (3.3) 

The only physical background to neutrinoless double (3 decay is from double (3 decay with 
two neutrinos. 

Note the following points: 
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• Since neutrinoless double (3 decay is a AL = 2 process, it is sensitive only to Majorana 
neutrino masses. Dirac masses conserve lepton number and do not contribute to this 
decay. 

• The neutrinoless double (3 decay rate is related to the neutrino mass-squared. It is 
also proportional to some nuclear matrix elements. Those matrix elements introduce 
theoretical uncertainties in extracting the neutrino mass from the signal, or in deriving 
a bound if no signal is seen. 

• Neutrinoless double (3 decay is sensitive to any AL = 2 operator, not only to the 
neutrino Majorana masses. Thus, the relation between the Neutrinoless double (3 
decay rate and the neutrino mass is model dependent. 

• The best bound derived from neutrinoless double (3 decay is m v < 0.34 eV [9]. 
C. Neutrino vacuum oscillation 

We now turn to the derivation of the neutrino oscillation formalism. We use several 
assumptions in this derivation. We emphasize, however, that in all practical situations it 
is correct to use these assumptions and that more sophisticated derivations give the same 
results (see, for example, [10]). 

In an ideal neutrino oscillation experiment, a neutrino beam is generated with known 
flavor and energy spectrum. The flavor and the energy spectrum is then measured at some 
distance away. If the flavor composition of the beam changed during the propagation it is a 
signal of neutrino oscillation, which indicates neutrino masses and mixings. 

The flavor of the neutrino is identified via its charged current interaction. It is convenient 
to express a flavor state in terms of mass states 



where Greek indices runs over the flavor states (a = e, /i, r) and Roman ones over the mass 
states (i = 1,2,3). An initially produced flavor state \u a ) evolves in time according to 



where we assume that all the components in the neutrino wave packet have the same mo- 
menta (pi = p). For ultra relativistic neutrinos we can use the following approximation 



(3.4) 



a 



(*)> = £e-^h>, 



(3.5) 




(3.6) 
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Then, the probability to observe oscillation between flavor a and (3 is given by 

Paf3 = \{vp\v a (t))\ 2 = 5 a p - ^Y^UaiUpiUajUpj] Sin 2 Xij, (3.7) 

i<j 

where in the last step we assumed CP conservation and 

AmU , 
*V = (3-8) 

In many cases the three flavor oscillation is well approximated by considering only two 
flavors. Then, 

/ cos 6 sin 9 \ 
U=[ , (3.9) 

\ — sm cos ) 

and for a ^ f3 we obtain 

P a p = sin 2 20 sin 2 x . (3.10) 
Since the neutrinos are relativistic we usually use L = t, E = p and write 

such that L osc is the oscillation length. The following formula is also convenient 

/Am 2 \ / L \ /GeV\ 

'"^bJlbJlT)' (312) 

The oscillation master formula, eq. (3.10), teaches us the following: 

• Longer baseline, L, or smaller energy, E, is needed in order to probe smaller Am 2 . 

• When L <^ L osc we can approximate sin 2 x ~ x 2 . In that case the signal is usually too 
small to be detected. 

• When L 3> L osc the energy spread of the beam and decoherence effects cause the 
oscillations to average out. That is, in our formula we should take the average value: 
(sin 2 x) = 0.5. In that case only a lower bound on the mass-squared difference between 
the two neutrinos can be obtained. 



D. Matter effects 

When neutrinos travel through matter the oscillation formalism is modified. While neu- 
trinos can scatter off the medium, in almost all relevant cases the scattering cross sections 
are very small, and the effect of scattering is negligible. The effect we are concentrating 
on is that of the forward scattering of the neutrinos. Like photons, when neutrinos travel 
inside a medium they acquire effective masses. There is no energy or momentum exchange 
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between the neutrinos and the medium. The effect of the medium is that it induces effective 
masses for the neutrinos. 

First we discuss the case of a medium with constant density. Moreover, we consider 
only matter at low temperature, which implies that it consists of electrons, protons and 
neutrons. Note that there are no muons, taus and anti-leptons in that case. The charged 
current interaction between the electron neutrinos and the electrons in the medium induces 
an effective potential for the neutrinos [1 1] 

V c = V2G F N e » 7.6 Y e ( P ) eV, (3.13) 

ylO^g/cnr 5 / 

where N e (N p , N n ) is the number density of electrons (protons, neutrons) and Y e = N e / (N p + 
N n ) is the relative electron number density. To get a feeling for the size of this potential, 
note that at the Earth core p ~ 10 g/cm 3 , which gives rise to Vc ~ 10~ 13 eV. 

As we discuss below, the current data indicates that m v > 10~ 3 eV, and thus m v ^> Vc- 
This seems to suggest that matter effects are irrelevant. This is, however, not the case since 
the matter effects arise from vector interactions while masses are scalar operators. The right 
comparison to make is between m 2 u and EVc where E is the neutrino energy. Since E ^> m v 
matter effects are enhanced and can be important. 

To explain this enhancement, we consider a uniform, unpolarized medium at rest. (Dis- 
cussion about the general case which includes all types of interactions and arbitrary polar- 
ization can be found, for example, in [12].) In that case the neutrino feels the four-vector 
interaction V^ = (Vc, 0,0,0). Due to Vc the vacuum dispersion relation of the neutrino, 
VnP^ = m2 ' * s m °dified as follows 

2 

777, 

(p l *-Vj(p> t -V> t )=m 2 E^p + V c + — , (3.14) 

zp 

where the approximation holds for ultra relativistic neutrinos. We learn that the effective 
mass-squared in matter, m^, is given by 

m 2 m = m 2 + A, A = 2EV C , (3. 15) 

where we used p ~ E. It is the vector nature of the weak interaction that makes the matter 
effects practically relevant. 

Interactions that are flavor universal only shift the neutrino energy by a negligibly small 
amount and do not affect the oscillation. Non-universal interactions, however, are important 
since in their presence the values of the effective masses and mixing angles are different from 
their vacuum ones. While the weak neutral current is flavor universal, the charged current 
is not. In normal matter, where there are electrons but not muons and taus, only electron 
neutrinos interact via the charged current with the medium. The mixing matrix is modified 
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by this non-universal matter effect such that the effective squared mass difference and mixing 
angle are given by 



Am 2 m = ^(Am 2 cos2#- A) 2 + (Am 2 sm29) 2 , (3.16) 

Am 2 sin 26* 

tan2# m = — — — -, 

Am 2 cos 29 — A 

where the subindex m stand for matter. The oscillation probability is then obtained from 
(3.10) by replacing the vacuum masses and mixing angles by the corresponding parameters 
in matter 

Am 2 L 

P a/3 = sin 2 29 m sin 2 x m , x m = ™ . (3.17) 

2hi 

The following points are worth mentioning regarding eq. (3.16): 

• The vacuum result is reproduced for A = 0, as it should. 

• Vacuum mixing is needed in order to get mixing in matter. 

• To first order in x m the matter effects cancel in the oscillation probability. To see this 
note that x sin 29 = x m sin2# m + 0(x^ n ). Therefore, when approximating sinx m ~ x m 
eq. (3.17) reduces to the oscillation probability in vacuum, eq. (3.10). 

• For Am 2 cos 29 ^> \A\, the matter effect is a small perturbation to the vacuum result. 

• For Am 2 cos 29 <C \A\, the neutrino mass is a small perturbation to the matter effect. 
In that case the oscillations are highly suppressed since the effective mixing angle is 
very small. 

• For Am 2 cos 29 = A, the mixing is maximal, namely it is on resonance. 
E. Non-uniform density 

When the matter density is not constant there are further modifications to the oscillation 
formalism. Density variation results in changing the effective neutrino masses and their 
mixing angles. Then, the flavor composition of the neutrinos along their path is a function 
of the medium density profile. 

At any point r on the neutrino path we define the derivative of the density 

= dMr)__ 
dr 

For constant density, A' = 0, the flavor conversion probability is controlled by the effective 
masses and mixing angles. For varying density, A' ^ 0, there are extra parameters that 
affect the flavor conversion probability. Most important is the adiabatic parameter 

_ Am 2 sin 2 2fl A^ 
Q{r) ~ Ecos29 A'(rY [ j 
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In the adiabatic limit, Q >> 1, the density variation is slow. In this case the transition 
between effective mass eigenstates is highly suppressed, and the constant density formalism 
can be applied locally. In the non- adiabatic limit, Q < 1, the density variation is fast. Then, 
transition between effective mass eigenstates is possible, and the constant density formalism 
cannot be used. Both limits can be of interest in reality. 

The best known example of the effect of the density variation is the MSW effect [11, 13]. It 
can account for a very efficient flavor conversion in the case of small mixing angle. Consider 
a two generation model where m 2 > mi with small vacuum mixing angle, 9 <C 1. That is, in 
vacuum the muon neutrino is mainly z/ 2 and the electron neutrino is mainly v\. We consider 
neutrinos that are produced in the Sun core where the density is large. In particular, we 
are interested in the case where at the core of the Sun, r = 0, the matter induce effective 
electron neutrino mass-squared is much larger than the vacuum mass-squared difference, 
A(0) 3> Am 2 . Then, at r = flavor eigenstates are almost pure effective mass eigenstates 

m {r = 0)^tt/2, v e {r = 0) w z/ 2 m , v^r = 0) « <. (3.20) 

(By writing u a v™ we mean that a very small rotation in flavor space is needed in order 
to move from the flavor eigenstate v a to the effective mass eigenstate vf 1 .) In particular, the 
produced electron neutrino is almost a pure v 2 . Also at the edge of the Sun, r = Rq, the 
flavor eigenstates are almost pure mass eigenstates 

9 m (r = Rq) -> 9 < 1, u e (r = Rq) « u u v^r = Rq) « v 2 . (3.21) 

When the adiabatic limit applies, neutrinos propagate in the Sun as effective mass eigen- 
states. Since the neutrinos are produced as almost pure u 2 , they leave the Sun as v 2 . Since 
in vacuum v 2 is almost a pure z/ M , we learn that there is almost full conversion from v e to v^. 
This phenomenon, which is called MSW resonance conversion, is demonstrated in fig. 2. 

In the more general case, corrections to the adiabatic limit are taken into account. In 
particular, transitions between the effective mass eigenstates are important in the vicinity 
of the resonance. Such level crossing is described by the Landau-Zener probability, Plz- 
Then, the flavor conversion probability is given by the Parke formula [14] 

Pee = - [1 + (1 - 2P LZ ) cos 29 m cos 29] , (3.22) 

where 9 m is the effective mixing angle at the production point. In this general case, the 
flavor transition probabilities become a rather involved function of the neutrino and medium 
parameters. It is this richness that makes the flavor composition of the solar neutrinos a 
complicated and not a monotonic function of their energies. As it is shown below this is 
important in order to be able to explain the solar neutrino data. 
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FIG. 2: A demonstration of an MSW resonance conversion. 9 m is the effective mixing angle in 
matter, A is the effective mass due to the matter effect defined in (3.15) and Ar is the resonance 
point. The neutrinos are produced at a region where A > Ar. They propagate as effective mass 
eigenstates, which is shown in the figure as solid lines. In the plot the effect is demonstrated for 
two different values of the vacuum mixing angles. (The plot is taken from [4].) 

F. Neutrino oscillation experiments 

Ideally, neutrino oscillation experiments measure the neutrino flavor transition probabil- 
ities as a function of the neutrino energies and their path. In practice, however, there are 
many experimental complications. The basic setup of any neutrino oscillation experiment 
is as follows: Neutrinos are produced and then propagate until few of them are detected at 
another location, far away from their production point. Therefore, in order to be able to 
probe the fundamental physics parameters we would like to know and control the following: 

• The total flux, the flavor composition and the energy spectrum of the beam. 

• The distance traveled by the neutrinos and the matter profile that they passed. 

• The total flux, the flavor composition and the energy spectrum of the detected neu- 



In reality, it is impossible to know and control all of these parameters. Therefore, a lot of 
work is involved in order to get the maximum out of each experimental setup. 

In general we distinguish between two kinds of experiments. Disappearance experiments 
search for reduction in the flux of a specific neutrino flavor 



trinos. 



N[u a (L)} < N[u a (0)} 



(3.23) 
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where N[u a (L)] is the number of neutrinos of flavor a in the neutrino beam at distance L. 
This would imply 

Paa{L)<\. (3.24) 

Appearance experiments, on the other hand, look for enhancement of a flux of a specific 
flavor 

N[v p {L)\ > N[vp(0)} , (3.25) 

which would imply 

Pa(s(L) > 0. (3.26) 

In particular, in many cases iV[z/g(0)] = 0. Then, any observation of neutrinos with flavor f3 
is an appearance signal. 

Oscillation experiments are usually named after their production place and mechanism. 
In the following we discuss solar neutrinos, atmospheric neutrinos and terrestrial neutrinos. 

IV. SOLAR NEUTRINOS 

Solar neutrinos receive a lot of attention since there are indications for disappearance 
of electron neutrinos. These indications are collectively referred to as the solar neutrino 
problem or the solar neutrino anomaly. The best explanation for the solar neutrino problem 
is neutrino flavor oscillation. 

The use of the words "problem" or "anomaly" should not worry us. They are used only 
to indicate the fact that solar neutrinos pose a problem to the SM. We think, however, 
that the SM is an effective theory of Nature, and thus, the fact that it faces problems is 
not surprising. Actually, we expect that neutrinos have masses and that they mix, and 
consequently, that the solar neutrino data cannot be explained by the SM. 

A. Solar neutrinos production 

Based on our understanding of stellar evolution we know that nuclear reactions fuel the 
Sun. Nuclear processes involve only the first generation of fermions and thus only electron 
neutrinos are generated in the Sun. Moreover, since the Sun is burn by fusion reactions it 
produces only neutrinos and not anti-neutrinos. The neutrino energies are relatively small, 
E v < 10 MeV. Therefore, only disappearance experiments can be done; the neutrinos 
energies are too small to produce muons or taus in any target at rest. 

There are many processes that produce solar neutrinos. The main reaction chain, which 
produces about 98.5% of the solar energy (and most of the solar neutrinos) is the pp cycle 

Ap -> 4 He + 2e + + 2z/ e + 2 7 . (4.1) 
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FIG. 3: The Bahcall and Pinsonneault 2000 solar model prediction of the neutrino flux as function 
of the neutrino energy. This figure was taken from [15] , where many more details on solar neutrinos 
can be found. 



The neutrinos emitted in this cycle have continuum spectrum which goes up to E v < 
0.42 MeV. Assuming that the Sun is in equilibrium, one can use its current luminosity 
to calculate the neutrino flux from the pp reaction chain. 

The other nuclear reactions in the Sun are not very important for its luminosity, but they 
are important for neutrino physics. The difference originates from the fact that photons are 
thermalized almost immediately after they are produced. Neutrinos, on the other hand, do 
not interact in the Sun, and leave the Sun with their original energies. Several sub-dominant 
reactions produce neutrinos with energies above the threshold of the pp chain. Since it is 
easier to detect high energy neutrinos, in many cases the experimental signals consist only 
of neutrinos that were produced by sub-dominant reactions. 

While the spectra of the different reactions are known, the sub-dominant reaction chains 
fluxes are not known in a model independent way. Solar models, namely, theoretical ar- 
guments, are used to obtain the total neutrino spectrum. As an example, the Bahcall and 
Pinsonneault 2000 (BPOO) solar model prediction for the neutrino flux is given in fig. 3. 
Since solar models involve theoretical inputs, one would like to perform experiments where 
the interpretation of the data can be done with minimum dependence on a specific solar 
model. 
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B. Solar neutrinos propagation 



Solar neutrinos are produced in the vicinity of the solar core. Then, they freely propagate 
out of the Sun and some of them travel to our detectors. During the night, neutrinos reach 
the detector after they have passed through the Earth. It is convenient to separate the solar 
neutrino propagation into three regions: inside the Sun, between the Sun and the Earth, 
and inside the Earth. We discuss each of these below. 

As already mentioned before, matter effects in the Sun can be important for flavor con- 
version. To calculate these effects we need to know the matter profile in the Sun and the 
rates of the different reactions in the different parts of the Sun. While these parameters are 
roughly known, there are uncertainties in their values. Therefore, we have to count on solar 
models to provide the needed inputs. 

The distance to the Sun is known to high precision. It varies by about 7% during the 
year. If vacuum oscillation is relevant to solar neutrinos then this distance variation can be 
used in order to get better control on the neutrino parameters. 

During the night neutrinos pass through the Earth on their way to the detector, while 
during the day they do not. Therefore, the Earth matter effect, which shows up as a day- 
night difference of the detected neutrino flux, provides further sensitivity to the neutrino 
parameters. Moreover, depending on the location of the detector and the time of the year 
the neutrino may pass through both the Earth mantle and its core or only through the 
mantle. This is important since the matter density in these two regions is different, and 
therefore can affect the oscillations [16]. 

C. Solar neutrinos detections 

There are three reactions that are used to detect solar neutrinos. The Charged Current 
(CC) interaction 

v e + n -> p + e~ , (4.2) 

where n is a neutron, can be used to detect only electron neutrinos. All the neutrinos can 
undergo Elastic Scattering (ES) 

vi + e~ — > v t + e~ . (4.3) 

This reaction is mediated by neutral current Z exchange for all neutrino flavors. For electron 
neutrinos there is also contribution from a W exchange amplitude. Due to this difference, 
the elastic scattering cross section for electron neutrinos is about six times larger than for 
muon and tau neutrinos. Finally, the cross section for the Neutral Current (NC) interaction 

Vi + n — > vt + n , (4.4) 
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is the same for all neutrino flavors (here n is a neutron). 

When presenting results of neutrino oscillation experiment it is convenient to define 

R = N obs /N MC , (4.5) 

such that N b s is the number of detected neutrinos and N MC is the number of events predicted 
from a solar model based Monte Carlo (MC) simulation assuming no flavor oscillation. For 
detectors that use neutral current reactions R should be one regardless if there are active 
neutrino flavor oscillations or not. For charged current and elastic scattering based detectors, 
however, R — 1 is expected only if there are no neutrino flavor oscillations, and R < 1 is 
expected if there are oscillations. 

There are several solar neutrino experiments that use the above mentioned basic reactions. 
They have different experimental setups that us various reactions with different sensitivities 
and thresholds. Next we describe these experiments. 

D. Chlorine: Homestake 

The first detection of solar neutrinos was announced 35 years ago. The neutrinos were 
detected in the Homestake mine in South Dakota [17]. The experiment uses the 

v e + 37 CI -> 37 Ar + e" , (4.6) 

charged current interaction in order to capture the electron neutrinos. 
We note the following points: 

• The signal is extracted off-line. Namely, only integrated rates are measured. Therefore, 
no day-night and almost no spectral information are available. (The only information 
about the neutrino spectrum is that the energy of the neutrino is above the threshold.) 

• The energy threshold is 0.814 MeV. Therefore, the experiment is sensitive mainly to 
the 7 Be and 8 B neutrinos, pp neutrinos cannot be detected in this experiment. 

• After many years of data collection the published result is [18] 

R^0.3. (4.7) 

Of course, in order to perform a fit to the neutrino parameters the above crude ap- 
proximation should not be used. Yet, quoting the order of magnitude, as we did in eq. 
(4.7) is sufficient to understand the main ingredients of the solar neutrino problem. 

The fact that the Homestake experiment found R < 1 was the first indication for solar 
electron neutrino disappearance. This result gave the motivation to build different types of 
detectors in order to further study the solar neutrinos. 
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E. Gallium: SAGE, GALLEX and GNO 



The Homestake experiment had relatively high threshold. The SAGE, GALLEX [19] and 
GNO [20] Gallium experiments were built in order to detect neutrinos with much smaller 
energies. In particular, their threshold is low enough to detect pp neutrinos. In the Gallium 
experiments the neutrinos are detected using the 

z/ e + 71 Ga^ 71 Ge + e- (4.8) 

charged current reaction. 

We note the following points: 

• Like the Chlorine experiment, the signal is extracted off-line. Thus, also here no 
day-night and spectral information are available. 

• The energy threshold is 0.233 MeV. Thus, Gallium detectors are sensitive to pp neu- 
trinos. The 7 Be and 8 B neutrinos also contribute significantly to the neutrino signal. 
Therefore, the total predicted flux is not solar model independent. 

• Combining the results of the three experiments, one obtains 

R « 0.6. (4.9) 

We learn that like in the Chlorine experiment, also here, the number of observed electron 
neutrinos is less than predicted. Moreover, the measured R is different for the two type of 
experiments. Since these two types are sensitive to different energy ranges, this fact indicates 
that the electron neutrino flux reduction must be energy dependent. 

F. Water Cerenkov: Kamiokande and SuperKamiokande 

Both the Chlorine and Gallium experiments extract their signal off-line. The Kamiokande 
and SuperKamiokande [21] water Cerenkov detectors can extract their signal on-line. There- 
fore, in these experiments day-night and spectral information can be extracted. The price to 
pay is that the energy threshold is much higher compared to the thresholds of the Chlorine 
and Gallium experiments. 

In the water Cerenkov detectors the neutrinos are detected by elastic scattering of the 
neutrinos off the electrons in the water 

v a + e- -> v a + e-. (4.10) 

We note the following points: 
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• Since the detection is done on-line, the time of each event is known, and also some 
information about the direction and energy of the neutrino in each event can be ex- 
tracted. These features enable the experiments to clearly show that the detected 
neutrinos are coming from the Sun. 

• The water Cerenkov detectors threshold is about 6 MeV. Thus, the Kamiokande and 
SuperKamiokande experiments are sensitive mainly to 8 B neutrinos. 

• Since elastic scattering is used for the detection, all neutrino flavors contribute to 
the signal. This fact is important if there are active flavor oscillations. Then, the 
converted neutrinos also contribute to the signal. This is in contrast to detectors that 
use charged current interactions where only electron neutrinos generate the signal. 

• The integrated suppression is found to be 

i2w0.5. (4.11) 
Spectrum and day-night information were also reported. 
We learn again that there is an energy dependent suppression of the electron neutrino flux. 



G. Heavy Water: SNO 

We already mentioned that neutrino flavor conversion is the most plausible explanation to 
the observed energy dependent electron neutrino flux suppression. Thus, it is interesting to 
measure solar neutrinos using neutral current reactions since in that case all neutrino flavors 
contribute with the same strength. Namely, neutral current based experiments measure the 
solar neutrino flux independent of whether there are active flavor oscillations or not. 

The SNO detector [22] was built for this purpose. Its apparatus contains an inner tank 
filled with heavy water and an outer tank filled with regular water. Neutrino can be detected 
by all the three basic reaction types 

CC(E th = 2.23 MeV) : v e + 2 H -> p + p + e~ , 
ES(E th ~ 6 MeV) : v a + e~ -> v a + e~ , 

NC(E th ~ 6 MeV) : u a + 2 H -> n + p + u a , (4.12) 

where E th is the threshold energy. Note that the presence of Deuterium ( 2 H) in the heavy 
water is crucial for the charged current and neutral current measurements. 
Note the following points: 

• The signal is extracted on-line. Thus day-night and spectrum information is available. 
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FIG. 4: The SNO data is plotted as three independent measurements of the muon or tau neutrino 
flux {4>ut) as a function of the electron neutrino flux (c/> e ). The three bands represent the neutral 
current, charged current and elastic scattering based measurements. The dotted line is the BPOO 
solar model prediction for the neutrino flux. It is plotted only in conjunction with the neutral 
current measurement since only this measurement is insensitive to electron to muon or tau flavor 
conversion. (This plot is taken from [23].) 



We define the ratios 

Re^, R N = ^, (4.13) 

where Rcc (Res, Rnc) is the ratio between the observed number of charged current 
(elastic scattering, neutral current) events and the predicted number assuming no 
oscillation. Extracting Re and Rn has several advantages. First, many systematic 
uncertainties cancel in their measurements. Second, their predicted values are almost 
solar model independent and depend only on the neutrino parameters. In particular, 
Re 7^ 1 or R^ ^ 1 can occur only if electron neutrinos are converted to muon or tau 
neutrinos. 

SNO found Re ^ 1 and Rn ^ 1. This result confirms the existence of or v T 
component in the solar neutrino flux at 5.3a. In addition, the BPOO solar model 
prediction is confirmed by the neutral current measurement. The SNO result [23] is 
illustrated in fig. 4. 
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H. Solar neutrinos: fits 



The general picture emerging by combining all the solar neutrino data is as follows. All 
charged current and elastic scattering based measurements found less events than predicted 
for massless neutrinos. This suppression is energy dependent. The BPOO solar model predic- 
tion for the neutral current based measurement, which holds also if there are active neutrino 
flavor conversions, was confirmed by SNO. 

The first and most robust conclusion is that the SM, which predicts massless neutrinos, 
fails to accommodate the data. While there are several exotic explanations that may be 
able to explain the data [24], the simplest and most plausible one is that of active neutrino 
flavor oscillation 

v e ^>v x , x = (x,t. (4.14) 

In that case, experiments that use reactions that are flavor dependent, namely charged cur- 
rent interactions and elastic scattering, are expected to detect fewer events than anticipated. 
On the other hand, experiments that use neutral current interactions, which are flavor blind, 
should not observe such reduction. Indeed, this is what the data shows. 

Assuming that the neutrinos are massive, the data can be used to determine the neutrino 
masses and mixing angles. While three flavor oscillation fits are available [25] , it is a good 
approximation to fit the data assuming only two neutrino mixing. An example of such fit is 
shown in fig. 5. (Recent fits can be found in [25]. Note, however, that whenever new data 
is published, new fits are performed and published usually within days.) The regions where 
a good fit is found are indicated in the figure. The best fit is obtained in the Large Mixing 
Angle (LMA) region. The recent SNO data actually exclude the Small Mixing Angle (SMA) 
region. The LOW and the Vacuum Oscillation (VO) regions, are still allowed by all solar 
data. They are excluded, however, by the recent KamLAND result, see below. 

To conclude, there are strong indications for active v e — > v x oscillation. This conclusion 
is robust since it has practically no dependence on the solar model. The best fit to the 
neutrino parameters is in the LMA region with 

Am 2 = 6 x 1(T 5 eV 2 , tan 2 9 = 0.4. (4.15) 

Future data and new experiments will teach us more about solar neutrinos. Besides having 
more statistics, future experiments will measure observables where we have only very little 
information at present. In particular, we expect to have more information about the day- 
night effect and the low energy neutrino spectrum. 
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FIG. 5: An example of a neutrino parameters fit to some solar neutrino data. The names of the 
regions where good fit is found are indicated in the figure. Note that by now this fit is outdated, 
and it is shown here for illustration only. (This plot is based on a plot taken from [26].) 



V. ATMOSPHERIC NEUTRINOS 



In the reaction chain initiated by cosmic ray collisions off nuclei in the upper atmosphere 
neutrinos are produced. These neutrinos are searched for by atmospheric neutrinos exper- 
iments. There are indications for muon neutrinos disappearance. These indications are 
called the atmospheric neutrino problem or the atmospheric neutrino anomaly. The likely 
explanation of this phenomenon is active neutrino oscillation. 

A. Atmospheric neutrino production 

Cosmic ray collisions off nuclei in the upper atmosphere produce hadrons, mainly pions. 
In their decay chain 

7T+ -> yL + p^ fi + -> e + v e V^ (5.1) 
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(and the charge conjugated decay chain) neutrinos are produced. Note that almost exclu- 
sively electron and muon neutrinos are produced. Based on (5.1) we expect 



where N(v e ) is the number of neutrinos of flavor t. There are corrections to this prediction. 
Other mesons, mainly kaons, are also produced by cosmic rays. Their decay chains produce 
different neutrino flavor ratio. Another effect is that at high energies the muon lifetime in 
the lab frame is much longer than in its rest frame. In fact, it is long enough such that not 
all the muons decay before they arrive to the detector. All in all, the neutrino production 
rates and spectra are known to an accuracy of about 20%. The ratio N(is^)/N(is e ) is known 
better, to an accuracy of about 5%. 

The atmospheric neutrino energies are relatively large, E v > 100 MeV. Therefore, both 
disappearance and appearance experiments are possible. In particular, since the initial tau 
neutrino flux is tiny, searching for tau appearance is an attractive option. 

B. Atmospheric neutrino propagation 

Since neutrinos are produced isotropically in the upper atmosphere, they reach the de- 
tector from all directions. The neutrino zenith angle in the detector is correlated with the 
distance it traveled. Neutrinos with small zenith angle, namely those coming from above, 
travel short distances (of order 10 2 km). Neutrinos with large zenith angle, namely those 
coming from below, are neutrinos that cross the Earth and travel much longer distances 
(of order 10 4 km). Since the production is isotropic, the neutrino flux in the detector is 
expected to be isotropic in the absence of oscillation. (There are small known corrections 
due to geo-magnetic effects.) 

Matter effects can be important for neutrinos that cross the Earth. Oscillations that 
involve electron neutrinos or oscillations between active and sterile neutrinos are modified 
in the presence of matter. In contrast, the matter has almost no effect effect on z/ M — > v T 
oscillation. 

C. Atmospheric neutrino detection 

There are two types of atmospheric neutrino detectors. The water Cerenkov detectors, 
1MB, Kamiokande and SuperKamiokande [21], use large tanks of water as targets. The 
water is surrounded by photomultipliers which detect the Cerenkov light emitted by charged 
leptons that are produced in charged current interactions between the neutrinos and the 




(5.2) 
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FIG. 6: Summary of the experimental measurements of R denned in eq. (5.3). (This plot is taken 
from [4].) 

water. Compared to electrons, muons are likely to create sharper Cerenkov rings. Thus, the 
shape of the ring helps in determining the flavor of the incoming neutrino. Iron calorimeter 
detectors are based on a set of layers of iron which act as a target, and some tracing elements 
that are used to reconstruct muon tracks or showers produced by electrons. 

Both water Cerenkov and Iron calorimeter detectors identify the neutrino flavor via its 
charged current interactions. They have some ability to look for neutral current interactions. 
Information on the neutrino energy and direction can also be extracted by both methods. 



D. Atmospheric neutrino data 

All atmospheric neutrino experiments measure the double ratio 

N'i IN\ 

T3 OPS I obs (r o\ 

^MC/ ly MC 

where N l MC is the expected number of flavor I type events and N^ bs are the observed ones. 
The advantage of using this double ratio is that many systematic and theoretical errors 
cancel in this ratio. Almost all of the experiments found R < 1; see fig. 6. 

The observation of R < 1 can be explained by muon neutrino disappearance, electron 
neutrino appearance, or both. The SuperKamiokande experiment also measured the zenith 
angle dependence of the neutrino flux. The data indicate that the preferred explanation is 
that of muon neutrino disappearance. For example, fig. 7 shows that the low energy v e flux 
agrees while the flux disagrees with the MC prediction. 
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FIG. 7: The low energy electron and muon neutrino fluxes compared to the MC predictions as 
observed by SuperKamiokande. The solid lines are the expected fluxes assuming no oscillations. 
The dashed lines are the best fits to neutrino oscillation. The solid circles are the data points. 
(This plot is taken from [27].) 



E. The atmospheric neutrino problem 



The zenith angle dependence of the muon neutrino flux as well as R < 1 cannot be 
explained with SM massless neutrinos. This is called the atmospheric neutrino anomaly. 
In order to explain it disappearance of muon neutrinos is needed. While there are several 
exotic explanations [24] , the most attractive explanation is active neutrino flavor oscillation. 
In particular, the best fit is achieved for — > v T oscillation with 

Am 2 = 2.6 x 1(T 3 eV 2 , sin 2 IB = 0.97. (5.4) 

The hypothesis of neutrino oscillation can be further tested. First, terrestrial long base- 
line experiments probe the same region of parameter space as the atmospheric neutrino 
experiments. Indeed, as discussed in more detail below, the recent K2K long baseline exper- 
iment results agree with the atmospheric neutrino data. Second, search for tau appearance 
in the atmospheric neutrino data is also possible. The SuperKamiokande experiment found 
such indications [28], but they are not yet convincing. 



VI. TERRESTRIAL NEUTRINOS 

One disadvantage of solar and atmospheric neutrino measurements is that since the neu- 
trinos are not produced on Earth, there is no control over the neutrino source. Neutrinos 
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from accelerators and nuclear reactors have the advantage that their source is known. That 
is, we know to high accuracy the initial neutrino flux, the energy spectrum and the flavor 
composition of the neutrino source. Moreover, in some cases these parameters can be tuned 
in order to maximize the experimental sensitivity. 

The main advantage of solar and atmospheric neutrinos is that they are sensitive to very 
small Am 2 . That is, their sensitivity is much better compared to most terrestrial neutrino 
experiments. Therefore, it is not surprising that most terrestrial neutrino experiments did 
not find disappearance or appearance signals. Their data was used to set bounds on the 
neutrino parameters. There are, however, some exceptions. The sensitivity of the K2K [29] 
and KamLAND [30] experiments are similar to those of, respectively, the atmospheric and 
solar neutrino experiments. Recently, positive disappearance signals were found in these two 
terrestrial experiments. The situation is different with regard to the LSND [31] experiment. 
In that case there is a claim for an appearance signal. This signal can only be explained 
with neutrino masses that are much larger compared to the masses deduced from the solar 
and atmospheric neutrino data. Next we describe these three experiments and their results. 

A. K2K 

In order for accelerator neutrinos to be sensitive to the same range of parameter space 
as the atmospheric neutrino experiments, long baseline, of order 10 3 km, is needed. In such 
long baseline experiments it is possible to search for disappearance or v T appearance. 
In these experiments a neutrino beam from an accelerator is aimed at a detector located 
far away. The first operating long baseline experiment is the K2K experiment. It uses an 
almost pure beam that is generated at KEK and is detected at SuperKamiokande, which 
is about 250 km away. 

Recently, the K2K experiment announced their first result [32] . They observed 56 muon 
neutrino events with an expectation of about 80. The disappearance can be explained 
by Vfi — > v x flavor oscillation. The best fit is found for the following values 

Am 2 = 2.8 x 10" 3 eV 2 , sin 2 26 = 1.0. (6.1) 

These values of the neutrino parameters are very close to those indicated by the atmospheric 
neutrino data, see eq. (5.4). The statistical significance of the K2K oscillation signal is 
smaller than the atmospheric neutrino ones. The importance of the K2K result lies in the 
fact that it provides an independent test of the atmospheric neutrino parameter space. 
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B. KamLAND 



Reactor neutrinos have much smaller energies compared to neutrinos produced at ac- 
celerator. Thus, long baseline reactor experiments are sensitive to smaller values of Am 2 . 
Another important difference is that reactors generate only electron anti-neutrinos, while 
accelerators produce mainly muon neutrinos. Thus, reactor experiments can be used to 
provide independent probes of the solar neutrino parameter space. 

The KamLAND experiment is placed in the Kamioka mine in Japan. This site is located 
at an average distance of about 180 km from several large Japanese nuclear power stations. 
The initial fluxes and spectra of the V e emitted in each reactor are known to a good accuracy, 
since they are related to the power that is generated in each reactor. Thus, measurement 
of the total flux and energy spectrum of the V e at KamLAND can be used to search for 
disappearance. This is particularly interesting since the KamLAND setup is such that it is 
sensitive to the region of the parameter space indicated by the LMA solution to the solar 
neutrino anomaly. 

Recently, the KamLAND experiment announced their first result [33]. They found 

N h 

— — = 0.611 ±0.085 ±0.041 (6.2) 

Nmc 

where N f, s is the number of observed events and Nmc is the expected number. This result 
can be explained by V e — > V x neutrino flavor oscillation with the best fit at 

Am 2 = 6.9 x 10" 5 eV 2 , sin 2 29 = 1.0. (6.3) 

These values are very close to the ones indicated by the solar neutrino data, see eq. (4.15). 

The importance of this result is twofold. First, it provides a completely solar model inde- 
pendent test for the neutrino oscillation solutions of the solar neutrino problem. Moreover, it 
discriminates between the different possible solutions, pointing at the LMA as the favorable 
one. 



C. LSND 

The set up of the LSND experiment [31] is as follows. Neutrinos are produced by sending 
protons on a fixed target that generates pions. Almost all the negatively charged pions are 
absorbed in the target before they decay. The neutrinos are produced via 

7T+ -> y+v^ /i + -> e + v e V fJL . (6.4) 

Then, the neutrinos travel for about 30 meter to the detector. Since the beam contains u^, 
u e and V^, with negligible fraction of V e , it best to search for V e appearance. 
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The LSND collaboration announced a 3.3a indication for v e appearance 

P(V^ -> V e ) = (2.64 ± 0.67 ± 0.45) x 10~ 3 . (6.5) 

This result can be explained with z/ M — > V e oscillation where the best fit is achieved for 

Am 2 = 1.2 eV 2 , sin 2 2# = 0.003. (6.6) 

Note that the value of Am 2 in (6.6) is much larger than the values indicated by the solar 
and atmospheric neutrino data. 

The LSND indications for neutrino masses are not as strong as the solar and atmo- 
spheric neutrino ones. First, the statistical significance of the LSND signal is rather low 
(traditionally, only a 5a effect is called a discovery). Second, the LSND signal still needs 
an independent confirmation. The Karmen experiment [34], which has a setting similar to 
LSND, but with a shorter baseline, L ss 17 meter, did not find an appearance signal. Soon, 
the MiniBooNE experiment [35] will be able to clarify the situation. 

VII. THEORETICAL IMPLICATIONS 

As we have seen, there are solid experimental indications for neutrino masses. The 
most important implication of these results is also the most simple one: the SM is not the 
complete picture of Nature. Namely, there are experimental indications that the SM is only 
an effective low energy description of Nature. 

There are basically two ways to extend the SM. One way is to assume that there is new 
physics only at a very high scale, much above the weak scale. Then, neutrino masses can be 
accounted for by the seesaw mechanism. Alternatively, one can think about weak scale new 
physics. Neutrino masses that are generated by such new physics are, in general, too large 
and some mechanism is required in order to explain the smallness of the neutrino masses. 
In the following we concentrate on the former option. 

A. Two neutrino mixing 

We start by considering only one type of results at a time. Namely, we draw conclusions 
from the solar and KamLAND data or from the atmospheric and K2K data. In that case 
there is one theoretical challenge imposed by the data: how to generate neutrino masses at 
the right scale? 

The atmospheric neutrinos and K2K results indicate — > v T oscillation with 

Am 2 AN ~ few x 10" 3 eV 2 . (7.1) 
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The solar neutrinos and KamLAND results indicate v e — > v^ T oscillation with 

Am^ N ~ l(T 4 eV 2 . (7.2) 

We use (2.9) in order to obtain the scale of the new physics that is responsible for neutrino 

masses 

m„ = -£ =^ M = — — . 7.3 

M m v v ' 

In order to find the high energy scale, M, we need to know m v and mo- The most likely 
situation is that the neutrinos are not degenerate 

A < -1. (7 .4) 



rrif + m 2 



Then, we use m v ~ V Am 2 . In general, the Dirac masses are of the order of the weak 
scale times dimensionless couplings. The problem is that we do not know the value of these 
couplings. In the following we make two plausible choices for their values. Motivated by 
5*0(10) GUT theories, we assume that these couplings are of order of the up type quark 
Yukawa couplings. In particular, for the heaviest neutrino, the coupling is of the order of 
the top Yukawa coupling. Another plausible choice is to use the charged lepton Yukawa 
couplings as a guide to the neutrino ones. 

Considering atmospheric neutrino and K2K data we use m v ~ few x 10~ 2 eV and get 

m D ~ m t M ~ 10 16 GeV, 

m D ~m T =>• M ~ 10 11 GeV. 

If instead we use m u ~ 10~ 2 eV as indicated by solar neutrinos and KamLAND data, we 
obtain 

m D ~m t =^ M ~ few x 10 16 GeV, 
m D ~ m r =>. M ~ 10 12 GeV. 
We emphasize the following points: 

• Both eqs. (7.5) and (7.6) imply that there exists a new scale between the weak scale 
and M P1 ~ 10 19 GeV. 

• The values obtained assuming uid ~ m t in both eqs. (7.5) and (7.6) are very close to 
the GUT scale, M G ut ~ 3 x 10 16 GeV. Therefore, we can say that neutrino masses 
are indications in favor of unification. 

• In several new physics models, a new scale, usually called the intermediate scale, 
Mi n t ~ ^rnivMpi ~ 10 11 GeV, is introduced. For example, supersymmetry breaking 
has to occur at this scale if it is mediated via Planck scale physics to the observable 
sector. The values obtained assuming ~ m T in both eqs. (7.5) and (7.6) are very 
close to M[ n t. This could be an indication that neutrino masses are also generated by 
such models. 
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B. Three neutrino mixing 



When we combine the solar, KamLAND, atmospheric and K2K data, we have to consider 
the full three flavor mixing. There are three mass differences that control the oscillations. 
They are subject to one constraint 

Am 2 12 + Am 2 23 + Am 2 31 = . (7.7) 

The data indicates that the mass-squared differences are hierarchical, |Am 2 2 | <C |Am| 3 |. 
This implies that |Am| 3 | ~ jAm^], and therefore, that there are only two different oscilla- 
tion periods. 

There are three mixing angles. Solar neutrino oscillations depend mainly on 612 and 
atmospheric neutrino oscillations depend mainly on # 23 . The third angle, #13, is known 
to be small, mainly from terrestrial experiments (see, for example, [4] for details). Thus, 
the oscillation phenomena are described by two mass-squared differences and three mixing 
angles 

Am 2 12 ~ 10" 4 eV 2 , Am 2 3 ~ few x l(T 3 eV 2 , (7.8) 
012 ~ 1, O23 ~ 1, 13 < 0.2. 

As before, we assume that the neutrinos are not degenerate. Then, we learn from (7.8) that 
the neutrino masses are somewhat hierarchical with two large and one small mixing angles. 

The theoretical challenge imposed by considering three neutrino flavor mixing is to explain 
the flavor structure of (7.8). There are basically two approaches to do so. One is called 
neutrino anarchy [36]. It assumes that there are no parametrically small numbers. The two 
apparently small numbers in the neutrino sector 

*•-§• ™ 

are assumed to be accidentally small. That is, numbers of order five are considered natural. 
Therefore, this mechanism predicts that #13 is close to its current upper bound. Conse- 
quently, it also predicts that CP violating observables, which are proportional to the product 
of all the mixing angles and the neutrino mass-squared differences, should also be close to 
their current upper bounds. 

The other option to explain (7.8) is to assume that there is an underlying broken flavor 
symmetry that controls the neutrino sector parameters. Such a symmetry is often assumed 
in order to explain the quark sector flavor parameters. Within this framework it is assumed 
that there is one small parameter, e, such that all the small parameters of the theory are 
functions of it. In particular, 

— ~e n «l, 0i 3 ~e m «l, (7.10) 
m 3 
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where n and m are some positive integers. In that case we often refer to observables like 
m 2 /m 3 and 13 as parametrically small numbers. 

While flavor symmetries can explain the observed hierarchies in the quark sector, it is 
not straightforward to extend the treatment to the neutrinos. The main reason is that in 
general large mixing angles come without mass hierarchies and vice versa. Explicitly, we 
consider the two generation case. The neutrino mass matrix is 




(7.11) 



and we distinguish the following three cases 

• a^b,c =>- sin 2 #^Cl, mi/m 2 < 1. 

• a, b, c ~ 1, det(m) ~ 1 ==>- sin 2 6 ~ 1, m 1 /m 2 ~ 1. 

• a,b, c ~ 1, det(m) <C 1 =>- sin 2 ^^!, m 1 /m 2 < ^l. 

We learn that in order to have large mixing with mass hierarchy the determinant of the 
mass matrix must be much smaller than the typical value of its entries (to the appropriate 
power). This is not a common situation in flavor models. In particular, this is not the case 
in the quark sector. 



C. Four neutrino mixing 



While the LSND data cannot be considered as a solid indication for neutrino masses, one 
may like to try to accommodate it as well. In that case the situation become much more 
complicated. The reason is that three different mass-squared differences are needed 

Am 2 SN ~ l(T 4 eV 2 , Am\ N ~ few x 10 _3 eV 2 , Am£ SND ~ 1 eV 2 . (7.12) 

With three neutrinos, however, the constraint (7.7) implies that there are at most two 
different scales. Therefore, at least four neutrinos are needed in order to accommodate all of 
these results. Experimentally we know that there are only three active neutrinos, and thus 
we need at least one additional light sterile neutrino field. 

There are two theoretical challenges: First, a sterile neutrino, by definition, is a fermion 
that is a singlet under the SM gauge group. As such, it can acquire very large mass since 
there is no gauge symmetry that forbids it. Thus, a sterile neutrino naturally has very large 
mass. There are, however, several ideas that produce naturally light sterile states [37]. 

Even if there are light sterile states, complicated patterns of masses and mixing angles 
are required in order to accommodate all the data. In fact, even adding a sterile state does 
not help in obtaining a good fit to the data [38]. In the following we do not address the four 
neutrino mixing any further. 
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VIII. MODELS FOR NEUTRINO MASSES 



As we saw there are three different kinds of theoretical challenges imposed by the neutrino 
data. Considering separately the solar and the atmospheric neutrino data, the challenge is 
to generate neutrino masses at the right scale. Combining them, we would like to find ways 
to generate the flavor structure, in particular, to explain mass hierarchy with large mixing. 
In order to accommodate also the LSND result, we should find a mechanism that generates 
light sterile neutrinos, and construct even more complicated flavor structure. 

The simplest way to generate neutrino masses is the seesaw mechanism. It is particularly 
attractive when it appears in models that are well motivated because they address various 
theoretical puzzles. Below we describe two models where neutrino masses are generated as 
one aspect of a model. (See [2, 4] for more examples.) 

A. Grand Unified Theories (GUTs) 

There are several supporting evidences for supersymmetric grand unification. First, in 
the minimal supersymmetric SM (MSSM) the three gauge couplings unify at one point at a 
scale of 



Second, one generation of SM fermions fits into two SU (5) multiplets and into one multiplet 
for higher rank GUT groups. Neutrino physics provide further support for GUT, particularly, 
for GUTs that are based on SO (10) (or higher rank) gauge group. Next we describe two of 
the facts that make SO (10) GUTs attractive from the point of view of neutrino physics. 

The 15 degrees of freedom of one SM fermion generation fall into a 16 of SO (10). The 
one extra degree of freedom is a SM singlet. Thus, it can serve as a right handed neutrino 
field. Once SO (10) is broken, this singlet becomes massive. Consequently, the SM active 
neutrinos acquire small masses via the seesaw mechanism. Namely, neutrino masses are 
predictions of SO (10) GUT theories. 

SO (10) also predicts the rough scale of the neutrino masses. To see it recall that SO (10) 
relates the up type quark masses to the Dirac masses of the neutrinos. In particular, the 
Dirac mass of the heaviest neutrino is of the order of the mass of the top quark. The 
Majorana mass of the singlet is expected to be of the order of the GUT breaking scale, 
Mn = A Aqut where A is an unknown dimensionless number. Then, the seesaw generated 
mass of the heaviest active neutrino is of the order of 



3 x 10 16 GeV. 



(8.1) 



1713 ~ 



m\ 10 _3 eV 



(8.2) 



M N A 
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For A ~ O(10 _1 ) it is the mass scale that was found by solar neutrinos. Somewhat smaller 
A is needed in order to get the scale that is relevant for atmospheric neutrino oscillation. 
Generally, it is fair to say that SO (10) based theories predict neutrino masses in accordance 
with the experimental findings. 

B. Flavor physics 

The most puzzling features of the SM fermion parameters are their smallness and the 
hierarchies between them. These features suggest that there is a more fundamental theory 
where the hierarchies are generated in a natural way. Generally, a broken horizontal symme- 
try is invoked. (By horizontal symmetry we refer to a symmetry that distinguishes between 
the three SM generations.) 

We demonstrate this broken flavor symmetry idea by working in the framework of a 
super symmetric U(1)h horizontal symmetry [39]. We assume that the low energy spectrum 
consists of the fields of the MSSM. Each of the supermultiplets carries a charge under U(1) H . 
The horizontal symmetry is explicitly broken by a small parameter A to which we attribute 
a charge —1. Then, the following selection rules apply: Terms in the superpotential that 
carry a charge n > under H are suppressed by 0(A n ), while those with n < are forbidden 
due to the holomorphy of the superpotential. 

Explicitly, the lepton parameters arise from the Yukawa terms 

Z ■ 

YijLiEjH d + -j^-LiLjH u H u . (8.3) 

Here Lj (Ei) are the lepton doublet (singlet) superfields, H u (H d ) is the up type (down type) 
Higgs superfields, is a generic complex dimensionless 3x3 matrix that gives masses to the 
charged leptons, Z^ is a symmetric complex dimensionless 3x3 matrix that gives Majorana 
masses to the neutrinos and M is a high energy scale. The selection rule implies 

if (Li) + H(E j ) + H(H d ) =n =}► = 0(\ n ) 

H(L i ) + H(L j ) + 2H(H u ) = m =^ Z tj = 0(\ m ) (8.4) 

where we assume that the horizontal charges of all the superfields are non-negative integers. 
We learn that the dimensionless couplings that are naively of order unity are suppressed by 
powers of a small parameter. 

Next we demonstrate a realization of the neutrino anarchy scenario using the above 
framework. Consider the following set of U(1)h charge assignments 

H(H u ) = 0, H(H d )=0, (8.5) 
H(L ± ) = 0, H(L 2 ) = 0, H(L 3 ) = 0, 
H(E 1 ) = 8, H(E 2 ) = 5, H(E 3 ) = 3. 
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Then, the lepton mass matrices have the following forms 

/A 8 A 5 A 3 \ 



M l ~ (H d ) 



A 8 A 5 A 3 
\A 8 A 5 A 3 / 



M 



n 1 i\ 
1 1 i 

\i i i/ 



(8.6) 



We emphasize that the sign "~" implies that we only give the order of magnitude of the 
various entries; there is an unknown (complex) coefficient of 0(1) in each entry that we do 
not write explicitly. Eq. (8.6) predicts large mixing angles 



sin# 



23 



non-hierarchical neutrino masses 



sin6> 



12 



sin# 



13 



M ' 



m 9 ~ 



M ' 



and hierarchical charged lepton masses 

m e ~ (H d )\ s , m, ~ (H d )\ 5 , 



3 M ' 



m r ~ (iJ d )A 3 . 



•7) 



(8.8) 



(8.9) 



We see that this model predicts neutrino masses and mixing angles that are in accordance 
with the neutrino anarchy scenario. When we assume that both (H d ) and (H u ) are at the 
weak scale and that A ~ 0.2, this model also reproduces the order of magnitude of the 
observed charged lepton masses. 

Other, more complicated, flavor structure can also be generated using a similar frame- 
work. See [2, 4] for more details. 



IX. CONCLUSIONS 



Neutrino physics is important since it is a tool for probing unknown physics at very short 
distances. Since we know that the SM has to be extended, it is of great interest to search 
for such new physics. In general, new physics at high scale predicts massive neutrinos. 
Therefore, there is a strong theoretical motivation to look for neutrino masses. 

In recent years various neutrino oscillation experiments found strong evidences for neu- 
trino masses and mixing 

• Atmospheric neutrinos show deviation from the expected ratio between the fluxes of 
muon neutrinos and electron neutrinos. Moreover, the muon neutrino flux has strong 
zenith angle dependence. The simplest interpretation of these results is that there are 
z/„ — v T oscillations. 
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• The electron neutrino flux from the Sun is smaller than the theoretical expectation. 
Furthermore, the suppression is energy dependent. The total neutrino flux, as mea- 
sured using neutral current reaction by the SNO experiment, agrees with the theoret- 
ical expectation. The simplest interpretation of these results is that there are v e — u x 
oscillations, where v x can be any combination of and v T . 

• The K2K long baseline experiments found indications for muon neutrino disappear- 
ance. Moreover, the same neutrino parameters that account for the atmospheric neu- 
trino disappearance also explain the K2K data. 

• The KamLAND experiment found evidences for V e disappearance. The data provide 
an independent test of the same parameter space that is probed by solar neutrinos. 

The theoretical expectation and the experimental data fit comfortably. Yet, there are 
many unsolved problems associated with neutrino physics. In the future we expect to have 
more data and thus we will be able to learn more about neutrinos, and eventually, about 
the short distance new physics. 
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